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Abstract 

We show that the Riccati-Pade method is suitable for the calculation of the complex 
eigenvalues of the Schrodinger equation with a repulsive exponential potential. The 
accuracy of the results is remarkable for realistic potential parameters. 



1 Introduction 

Recently, Rakityansky et al [1] developed a method for the calculation of 
bound states and resonances of quantum-mechanical problems. It is based 
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on the approximation of the S-matrix by means of Pade approximants and 
the location of their poles. The authors mention that the methods for the 
calculation of bound states and resonances are usually developed separately 
in spite of the fact that those quantities share the same mathematical nature. 

Some time ago we developed the Riccati-Pade method that applies to bound 
states and resonances of separable quantum-mechanical models [2-11]. Al- 
though the RPM is less general than the approach proposed by Ratikyansky 
et al [1] it is nevertheless an interesting approach for comparison purposes and 
benchmark. 

The RPM [2-11] is based on a rational approximation to a modified logarith- 
mic derivative of the wavefunction that satisfies a Riccati equation. From the 
coefficients of the Taylor expansion of this logarithmic derivative, which are 
functions of the energy, we construct Hankel determinants. Their roots give 
rise to sequences that converge towards the bound states and resonances of 
the quantum-mechanical model as the determinant dimension increases. In 
most cases the rate of convergence is so great that the RPM yields extremely 
accurate real and complex eigenvalues. 

The rational approximation and the Hankel quantization condition appear 
to select square integrable functions and incoming or outgoing waves. Both, 
bound states and resonances emerge from the Hankel sequences because one 
does not introduce the boundary conditions at infinity explicitly. 

Earlier results for exponential potentials suggest that the RPM may not 
yield virtual states and that the Hankel sequences converge to a wrong limit, 
although, suspiciously close to the right answer, in the case of some reso- 
nances [9] . The purpose of this paper is to investigate this feature of the RPM 
more closely. 

In Section [2] we discuss a simple model that enables us to calculate the poles 
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of the scattering amplitude from the roots of Bessel functions. In Section [3] 
we apply the RPM to this model and compare approximate and exact results. 
Finally, in Section H] we summarize our results and draw conclusions. 



2 Model 

In this paper we test the performance of the RPM on the Schrodinger equation 

ij"{x) + - Ae-"^) ^{x) = (1) 

with the boundary condition ip{0) = 0. This model has proved useful in the 
past for the study of bound, resonance, and virtual states [12-14]. Besides, 
the exponential potential is a suitable representation of repulsive molecular 
interactions [13]. The change of variables q = ax, ^(g) = ip^q/a) leads to an 
eigenvalue equation with just one potential parameter: 

+ (e - Ae-") <l>(g) = (2) 
where e = E/a^ and A = A/a'^. 

A further change of variables z = 2V— Ae~''/^, Y{z) = $(2 ln[2V— A/2;]) trans- 
forms the eigenvalue equation ([2]) into the Bessel equation 



z' 



Y"{z) + zY\z) + - Y{z) = (3) 



where z/^ = — 4e. The general solution that satisfies the boundary condition 
at g = is [12] 



Y(z) = C 



J_. 2V-A JJz) - J, 2V-A JJz 



(4) 



If we assume that the RPM will provide the eigenvalues of those solutions 
that satisfy hmg^oo $(g) = 0, then it follows from the behaviour of the Bessel 
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function at origin Ju{z) ~ /^{^ + 1) that they should be roots of 



J, (2v^) = 



(5) 



with Re(z/) > 0. The roots of this equation are poles of the scattering ampli- 
tude [13]. 

The modified logarithmic derivative of an eigenfunction of equation ([21) 

fiQ) = - - ^ (6) 



satisfies the Riccati equation 



f{q) + ^ - f(qy + Ae-' -6 = (7) 



The Taylor series about the origin 



/(?) = E/.V\ /o = o 

i=o 



converges in a neighbourhood of g = and the coefficients fj depend on e. 

The main assumption of the RPM is that the roots of the Hankel determinants 
Hf){e) = 0, with matrix elements fj+j+d~i, i,j = 1, 2, . . . , -D, are suitable ap- 
proximations to the eigenvalues of the Schrodinger equation ([2]) [2-11]. Here, 
D = 2, 3, ... is the determinant dimension and d = 0,1, . . .. More precisely, 
we expect that there exists a sequence of roots e'^''^' of the Hankel determi- 
nants that converges to a given eigenvalue of that Schrodinger equation as D 
increases. 
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3 Results and discussion 



Previous applications of the RPM showed that the rate of convergence of the 
Hankel sequences is remarkable for both real and complex eigenvalues [2-11]. 
However, in the case of the exponential potential ([1]) the Hankel sequences were 
found to converge to a result slightly different from the one given by equation 
([5]) [9]. For example, Table[T]shows Hankel sequences for A = 0.5 and A = 2 and 
the corresponding exact results obtained from the quantization condition ([5]). 
Both Hankel sequences exhibit great convergence rate but they do not converge 
towards the expected result. We also appreciate that the disagreement between 
the exact and RPM eigenvalues increases as A decreases. In addition to the 
great convergence rate, the Hankel determinants exhibit clustering of roots 
about the limits of the sequences as D increases, which is an indication of 
satisfactory convergence and meaningful result. However, those limits do not 
completely agree with the exact results given by ([5]). 

Fig. [T] shows values of e(A) calculated by means of the RPM and equation 
([5]). As indicated above, the agreement between the RPM and exact results 
increases as A increases. For A < 0.4 the RPM yields complex eigenvalues in 
spite of the fact that the exact ones are real (see, for example, Atabek et al [13] 
for a more detailed discussion of this behaviour). Moreover, it is well known 
that the roots u of equation ([5]) tend to negative integers as A [13, 14] 
but we clearly see that the RPM eigenvalues do not exhibit this behaviour. 
However, it is most interesting that both real and imaginary parts of the 
RPM eigenvalues give a reasonable picture of the behaviour of the exact ones 
for large and moderate values of A as shown in Fig. [H In order to show the 
increasing agreement between the RPM and exact eigenvalue with A more 
clearly. Fig. [2] shows log — gR^A^j ^s a function of A. 

The discrepancy between the RPM and exact eigenvalues just mentioned is 
interesting from a mathematical point of view, but it is not a serious drawback 
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for physical purposes. It is well known that realistic potentials require much 
larger values of A than those in Fig. [1] [13]. If, for example we choose A = 100 
the rate of convergence of the Hankel sequence is much greater than the one 
in Table [1] and the limit agrees with the exact result to at least 20 digits (see 
Table [2]). Atabek et al [13] estimated the potential parameter for the 
repulsive state of BeH to be A ^ 134 that corresponds to A = 4489. In this 
case the rate of convergence of the Hankel sequence is even greater and we 
obtain the exact result to at least 20 digits with D = 15 and d = as shown 
in Table [2l It is worth mentioning that for such large values of A we find it 
easier to obtain the complex energies by means of the RPM than from the 
roots of the Bessel function. 

4 Conclusions 

We have shown that the Hankel sequences converge to a wrong limit that is 
close to the resonances of the repulsive exponential potential. At present we 
do not know the reason for this discrepancy that decreases as the potential 
strength increases. However the RPM eigenvalue as function of the potential 
parameter follows the trend of the exact one. Besides, for realistic potential 
functions that fit the interaction between molecular fragments [13], the RPM 
provides remarkably accurate results and it is probably more accurate than 
other methods. Unfortunately we are not aware of results for such great values 
of the potential strength, probably because other methods do not provide the 
complex eigenvalues so efficiently as the RPM. 

We should mention that the complex eigenvalues of the Schrodinger equation 
with the exponential potential ([T]) are not what physicist use to call resonances 
because the imaginary parts are too large. However, there has been some 
interest in their calculation with the purpose of reconstruction of the scattering 
amplitude from its poles [13, 14]. 
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Table 1 

Convergence of the Hankel sequences for two values of A 



A = 0.5 



D 


Re e 


Im e 


10 


-0.70545054582805260895 


0.26816598479688956569 


11 


-0.70545056661473611239 


0.26816596401669675062 


12 


-0.70545056816101389381 


0.26816596478162896896 


13 


-0.70545056810407142688 


0.26816596486012868610 


14 


-0.70545056805626480670 


0.26816596487152171555 


15 


-0.70545056805476309518 


0.26816596487187249577 


Ifi 


-0 7n54'in5fi8n54QQfi189^^ 


n 96816^964871 56659^41 


17 


-0 70545056805503003056 


26816596487158015859 


18 


-0.70545056805502881287 


0.26816596487157957123 


19 


-0.70545056805502843245 


0.26816596487157970586 


20 


-0.70545056805502836237 


0.26816596487157971713 


Exact 


-0.73985910415959609800 


0.24527511363052010569 


A = 2 


D 


Re e 


Im e 


10 


-0.66695560251365514674 


1.6211700286378446183 


11 


-0.66695559703717080232 


1.6211700273918467827 


12 


-0.66695559708065060107 


1.6211700264998970342 


13 


-0.66695559710083543597 


1.6211700264601278937 


14 


-0.66695559710850847182 


1.6211700264500121775 


15 


8 

-0.66695559710929174524 


1.6211700264509149572 



Table 2 

Convergence of the Hankel sequences for large values of A 



A = 100 



D 




Re e 




Im e 


5 


71, 


.535851840807002875 


37. 


,763655201763995538 


6 


71, 


.535265703486635320 


37. 


,763686464878442119 


7 


71, 


.535231112272632786 


37. 


,763673576025908553 


8 


71, 


.535229874257576952 


37. 


,763674377953883589 


9 


71, 


.535229860223483406 


37. 


,763674375421503193 


10 


71, 


.535229855328822688 


37. 


,763674374272030169 


11 


71, 


.535229855364554354 


37. 


,763674374318377973 


12 


71, 


.535229855364593976 


37. 


,763674374316154213 


13 


71, 


535229855364798246 


37. 


763674374316019149 


14 


71, 


.535229855364801073 


37. 


,763674374316016150 


15 


71, 


.535229855364801145 


37. 


,763674374316015925 


16 


71, 


.535229855364801148 


37. 


,763674374316015917 


17 


71, 


.535229855364801148 


37. 


,763674374316015916 



A = 4489 


D 


Re e 




Im € 



5 4158.9571348913180566 530.18487292989576140 

6 4158.9534557328433593 530.18603523832850612 

7 4158.9533491082024576 530.18608879147133066 

8 4158.9533468989068860 530.18609014361891762 

9 

9 4158.9533468479121394 530.18609015967783227 




Fig. 1. Eigenvalues of the exponential potential: Re e^^^ (empty circles), Im e^^^ 
(solid circles), Re e^^""^* (dashed line), Im e^^"'^* (solid line) 




Fig. 2. log I Re e^^""^* — Re e^^^\ (empty circles) and log |Im gexact e^^^\ (solid 
circles) 
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